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Vittorio Giovannetti, Diego Frustaglia, and Fabio Taddei 
NEST-CNR-INFM and Scuola Normale Superiore, 1-56126 Pisa, Italy 



00 

o 

O 
(N 

X) 
(D 

oo 



X3 

o 



> 
cn 

o 
o 



I 

o 
o 



X 
S3 



Rosario Fazio 

International School for Advanced Studies (SISSA), via Beirut 2-4, 1-34014 Trieste, Italy 
NEST-CNR-INFM and Scuola Normale Superiore, 1-56126 Pisa, Italy 
(Dated: February 18, 2008) 

We show that multi-mode entanglement of electrons in a mesoscopic conductor can be detected 
by a measurement of the zero- frequency current correlations in an electronic Hong-Ou-Mandel inter- 
ferometer. By this means, one can further establish a lower bound to the entanglement of formation 
of two-electron input states. Our results extend the work of Burkard and Loss [Phys. Rev. Lett. 
91, 087903 (2003)] to many channels and provide a way to test the existence of entangled states 
involving both orbital and spin degrees of freedom. 

PACS numbers: 03.67.-a,03.67.Mn,72.70.-f m,73.23.-b 



I. INTRODUCTION 



Differently from quantum optics, where entangled pho- 
tons are routinely detected by coincidence measurements, 
the issue of revealing the presence of entangled states in 
multi-terminal mesoscopic conductors (see, e.g., Ref. [l| 
and references therein) is still an experimental challenge. 
A number of theoretical proposal have been considered 
so far. Bell-like tests by means of measurements of 
current-fluctuations have been analyzed by several au- 
thors^. Quantum state tomography has been discussed 
in Ref. Alternatively, Burkard et al^ suggested that 
singlet and triplet electron pairs could give rise to devia- 
tions in the zero-frequency shot noise at the output ports 
of a 50/50 electronic beam splitter (BS) as compared to 
the value observed for an incoming beam of independent 
electrons. Following this work, the full counting statistics 
for entangled electrons in a BS has been studied in Ref. 
The effect of a Rashba spin-orbit term or a rotating mag- 
netic field in one of the incoming ports was discussed 
in Refs. and 0, respectively, and in Ref. H dephasing 
through additional reservoirs was included. While these 
works considered only singlet or triplet incoming states, 
in Ref. 9 the case of states generated in an Andreev dou- 
ble dot entangler at the input ports of the BS has been 
discussed. 

Our work is motivated by a recent paper by Burkard 
and Lossi^ where the authors derived a lower bound for 
the entanglement of formationii of arbitrary mixed spin 
states, through shot-noise measurements. In this paper 
we generalize their results to multi-mode input states by 
introducing the electronic analog of the Hong-Ou-Mandel 
(HOM)i^ optical interferometer. We assume that a two- 
electron (possibly entangled) input state is injected into 
a BS whose arms support many propagating modes, in- 
cluding both orbital and spin degrees of freedom. In this 
context, we first show that zero-frequency noise measure- 
ments allows one to detect entanglement for two-electron 
input states. Then, in the spirit of Ref. fl^ we derive a 
lower bound for the entanglement of formation in the case 
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FIG. 1: Schematics of the apparatus. For any energy E, two 
electrons enter from an entangler at ports 1 and 2. Each port 
consist of many channels, labeled by a spin (s) and orbital 
{I) index a = {s,£}. The BS (with transmissivity T that 
could be controlled by a gate voltage) mixes the incoming 
state, leading to an outgoing state propagating along ports 3 
and 4. Currents are measured at 3 and 4 for different values 
of a channel-dependent phase shift Lpa introduced in port 2. 
The presence of spin as well as orbital entanglement in the 
input state ([6]) is detected through the Fano factor ((TJ for 
current-noise correlations. 



of arbitrary number of incoming channels. This is rele- 
vant from an experimental point of view since it shows 
the possibility of characterizing entanglement in complex 
scenarios which go beyond the paragdimatic two-mode 
implementation of Burkard and Loss original scheme. 
Moreover, from a theoretical point of view, our analy- 
sis establishes a clear connection between the detection 
scheme of Ref. J0| and the HOM interferometer scheme 
proposed in optics. 

The main ingredient in the present approach is the 
study of the oscillations of the zero-frequency shot noise 
at the outgoing ports of a BS as a function of a con- 
trollable phase shift present in one of the incoming 
ports. As we shall see, this technique allows us to de- 
tect spin/orbital as well as more complex entangled input 
states. In view of its simple implementation we believe 
this scheme can offer a very convenient way for entangle- 
ment detection in mesoscopic multi-terminal conductors. 

The paper is organized as follows. In Sec. |lT] we in- 
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troduce the electronic analog of the HOM interferometer 
and, in Sec. lIIIi we demonstrate that it detects entangle- 
ment for two-electron input states. In Sec. llVl we discuss 
how the zero-frequency noise measurements can be re- 
lated to the entanglement of formation of the input state. 
The paper finally ends with conclusions and remarks in 
Sec. El 

II. THE HOM INTERFEROMETER 



In the HOM experimeniji^ pairs of multi-mode pho- 
tons propagating along two distinct optical paths inter- 
fere at a 50/50 BS after being phase shifted through con- 
trollable delays. The presence of quantum entanglement 
in the input state of the photons can be recovered by 
studying the coincidence counts at the output ports of 
the BS as a function of the relative delay between op- 
tical paths. In the electronic analog of the HOM inter- 
ferometer (c-HOM) metallic conductors play the role of 
photonic paths and current correlations that of coinci- 
dence counts. Apart from the expected differences re- 
lated to the statistics, the e-HOM interferometer main- 
tains the same entanglement-detection capability of its 
optical counterpart. 

The e-HOM is sketched in Fig.[TJ Pairs of electrons of a 
given energy E above the Fermi sea, prepared in a (possi- 
bly entangled) initial state, enter the interferometer from 
the input ports 1 and 2. Electrons passing through the 
port 2 undergo an additional, controllable, phase shift 
(white circle in Fig. [T]) before impinging on the BS. Zero 
frequency current correlations are measured at the out- 
put ports 3 and 4. The electron states, at energy E, are 
labeled by the indices (j, a) where j = 1, • • • ,4 labels the 
ports of the e-HOM interferometer, and where a is the 
composite index {£, s} with s referring to electron spin 
component along the quantization axis, and £ referring 
to the orbital channel. Following the Landauer-Biittiker 
scattering formulation of quantum transport ^'^ we intro- 
duce a set of fermionic operators for the incoming aj^a{E) 
and outgoing bj^a{E) states. They are connected via the 
scattering matrix by the relation 

b,AE)=Y,Sl'^,iE)a,,AE), (1) 
j' 

where we assume that the spin is conserved in the scat- 
tering process and that there is no channel mixing. The 
phase shifts at the input port 2 are described by the map- 
ping 

a2AE) ^ e'^'-a^AE) , (2) 

where (pa depends on some externally controlled param- 
eters. These transformations can be implemented by 
introducing local gate voltages and/or magnetic fields 
at the input port 2, allowing for an independent con- 
trol of orbital- and spin-related phases, respectively. 



In particular, by introducing local electric gatings Vi 
plus an effective magnetic anisotropy Bcff (due to, e.g., 
spin-orbit coupling®), for small perturbations we find 
= {eVe/Ei){keLo/2) + s{^iB^s/Ee){keLi/2), with h 
the electronic wave-number along the orbital channel £, 
El ~ (^^/2m)fc|, ^, the electronic magnetic moment, 
s = ±1 the spin, and Lq {Li) the length on which V( 
(Bes) acts. For simplicity we will assume a symmetric 
BS which does not suffer from backscattering. Hence the 
scattering matrix S has a block structure of the form 

s=[l 1)^s(''Ke), 

where «(")(£') describes the transformation of 
ai,a{E),a2.a{E) into &3, «(£'), 64,0 (-E) and it is ex- 
pressed by the matrix 

where T is the BS transmissivity that could be controlled 
by, e.g., a gate voltagci^ 

The current operator of the j-th port is defined asi^ 

E e-^-%jE)b,AE + ncu) 

E,u}.a. 

-a]jE)a,AE + na;)], (4) 

where v is the density of states of the leads and where 
a discrete spectrum has been considered to ensure a 
proper regularization of the current correlations. The 
zero- frequency current correlations are defined as 

- lim ^ r dh r dh {5I,ih)5IAt2)), (5) 

where the average (■ ■ • ) is taken over the incoming elec- 
tronic state, T is the measurement time and 6Ij = 
Ij — (Ij). We now show that by studying the functional 
dependence of Eq. ([5]) upon the phase shifts one can de- 
tect the presence of entanglement in the input state of 
the interferometer. 



III. ENTANGLEMENT DETECTION 

Consider the two- electron case in which, for a given 
energy E, one electron per port enters the interferometer 
from 1 and 2. The most general pure input state of this 
form can be expressed as 

I*) = n E 4jE)alAE) |0) , (6) 

E a, 13 

where the product is taken for energies in the range 
< E < eV, as though ports 1 and 2 were kept at a 
voltage bias V with respect to ports 3 and 4. Further- 
more, |0) is the Fermi sea at zero temperature and $q./5 is 
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the two-electron amplitude which we assume to be inde- 
pendent of E and satisfying the normalization condition 



l^a,/3| 



1. One can easily verify that, indepen- 



dently of $Q,/3, the average current is constant and equal 
to e^V/h. A straightforward calculation of the Fano fac- 
tor Fjji — Sjj' I (2e-\/ (/j) (/j')) leads to the expression 



F, 



34 



where 



/3,Q e 



(7) 



(8) 



is a (real) quantity which depends on the controllable set 
of phases {iy9a}. By using the Chauchy-Schwartz inequal- 
ity and the normalization condition of ^afi it follows that 
for generic input state ([6]) one has — 1 ^ u;$ ^ 1. How- 
ever, if \^) in Eq. ([S]) is a separable state with respect 
to the input ports 1 and 2, it is possible to show that 
is non negative. In fact, for a separable state ©, the 
two-electron amplitude factorizes as $a,/3.scp = Xol'x^'^^ 

with Xq ^ (x^^^ ) being the amplitude associated to the in- 
coming electron of port 1 (2). Replacing this expression 
in Eq. ([8]) we get 



Sep 



^ 



(9) 



By convexity the same result applies also to any separable 
mixed state of the form psep = Y.i Pi l*sGp(«))(*scp(i)l' 
with Pi > 0. Negative values of it;$ are hence a direct 
evidence of the presence of entanglement in the input 
state, i.eJ^ 



If ^^34<-T(l-r) 



The state is entangled. (10) 



Notice that the inverse implication does not hold: indeed 
there exist entangled states ([5]) which have positive 
and hence F34 ^ — r(l — T) (see Fig. [5]for an example). 
The quantity w$ therefore acts as an entanglement wit- 
nessi^ for the class of two-electron states analyzed here. 
The condition of Eq. pil)) has been derived in Ref. 10 for 
the spin entangled case. Here we extended its validity to 
a generic multi-mode entangled input state as defined in 
Eq. ([6]) . Equation (jlOp shows that a simple measurement 
of the Fano factor may be sufhcient to ascertain if a given 
state is entangled. Once the transmissivity T of the BS 
is known and the phases {fa} are tunable, the test pU]) 
is within current available experimental abilities. 

As an example we illustrate in Fig. [5] the detection of 
the entanglement between two-level electronic states with 
(pseudo)spin a at, e.g., the output of one of the 

spin or orbital entanglers so far proposed,- By defining 
$Q,,/3 = |^'Q,/3|e"'°'^ we rewrite the quantity in Eq. ^ as 



1$ 



T,Tl 



1$ 



(11) 



2|$|,x$x,||cos(A(^ + Ar;) 




FIG. 2: Entanglement detection scheme applied to the case 
of two- level electronic states with component a =Ti-l- The 
quantity w<s> is plotted as a function of Aip — tp^- ~ (p^^. A 
negative w<s> (shaded area) signals the entanglement in the 
input state ((6]). Results are shown for: (solid line) a detectable 
entangled state with $t,| = 1/^2, $x,t = ^i,! = 

^1,1 ~ l/\/T2i (dashed line) an undetectable entangled state 
with <1>|4 = $1,1 = $1,1 = ^3/10, <l>i,T = e~"''''VVTO; and 
(dotted line) a separable state with "l>u = — 1/2, and 



with Arj = rji^ — rjn and Aip = (p]- — <y3|. Entanglement 
can be guaranteed only for those input states leading to 
a negative w$ (shaded area in Fig. [2]) for some value 
of A(p. We show results for three typical states that can 
enter the BS analyzer: a detectable entangled state (solid 
line), an undetectable entangled state (dashed line), and 
a separable state (dotted line). 



IV. LOWER BOUND ON THE 
ENTANGLEMENT OF FORMATION 



Our analysis of the generic input state defined in Eq. 
proceeds further. By generalizing an argument presented 
in Ref. ilQ . in this section we show that the current noise 
measurement of Eq. ([7]) is related to the entanglement 
of formationii of the input state of the e-HOM inter- 
ferometer. Fundamental ingredients of our analysis are 
the d-dimensional generalized Werner state^^^ whose 
properties are briefly reviewed in the following for com- 
pleteness. 



A. Generalized Werner stated and generalized 
twirling transformations 

Consider a system composed of two d-dimensional sub- 
systems 1 and 2 and characterized by a separable canon- 
ical basis \al3) = \a)i (8) |/3)2 with a, /3 € {1, • • • , d}. Fol- 
lowing the notation of Ref. [l^ we introduce the general- 
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ized Werner statesi^ii^ 

2(1 - W) 



aw 



Q</3 



+ 



2W 



d{d-l) 



1/3 I 



, Q</3 



with We [0,1] and 



Q/3 



V2 



(12) 



(13) 



For d = 2 Eq. reduces to the standard definition 
of a qubit Werner statejiS The entanglement of forma- 
tionii Ef of the density matrix aw has been computed 
in Ref. |T^. Indeed one can verify that these states are en- 
tangled if and only if > 1/2 and that Ef{aw) = £iW) 
with £{x) being the function 



£{x) 



(here H{x) 




l+2y/x{l-x) 
2 



for X e [0,1/2] 
for X £ ]1/2,1] , 



(14) 



-xlOE 



(1 — x) log2(l — a;) is the bi- 



nary Shannon entropy). Another important feature of 
the Werner set (|12p is the fact that any input state p 
of the system can be transformed into one of the aw 
by means of Local Operations and Classical Communi- 
cations (LOCC), i.e. 



P 



(15) 



with W{p) being the expectation value of the observable 
p I on the input state p, i.e. 



w^(p) = E(*SH*i7>- 



(16) 



Q</3 



The map (|15l) is a generalized twirling transformation ^^ ' ' '^^ 
which can be expressed as follows, 



Twcr(p) = dU {U(E)U) p ([/■f (g) C/t) 



(17) 



with U being a generic unitary operator on the subsys- 
tem and dU being a proper measure in the space of such 
unitary transformations. It can be shown that Werner 
states aw are fix-points of Twor, i-e. T'wor(o'w) = o'w- 

By exploiting the properties of Twer, Eq. lfTS]) yields a 
lower bound to the entanglement of formation of the den- 
sity matrix p. In fact, since the entanglement of a state 
does not increase under LOCC transformations, we must 
have 



Efip) ^ Ef{awip->)=£{W{p)), 



(18) 



with £{x) as in Eq. Eq. (fT5|) shows that the quantity 
W{p) of Eq. can be used to bound the entanglement 



of p. In Sec. lIVBj we shall prove that the Fano factor ([7]) 
with all ^ provides a natural way of measuring the 
quantity (|16|1 for two-electron input states. 

The bound (|18p can be strengthen by noticing that 
the above derivation holds for all choices of the sepa- 
rable canonical basis |a/?). Indeed for any inequivalcnt 
definition of such basis one gets a new set of Werner den- 
sity matrices and a corresponding new LOCC twirling 
transformation. Therefore, since £{x) is non decreasing, 
Eq. (fT^ can be replaced by 



Efip) ^ £{W{p)) 



(19) 



where W{p) is the maximum value of the quantity in 
Eq. (|16p with respect to any possible choice of the sep- 
arable basis I a/3). To analyze the performance of the 
e-HOM set-up we can focus only on those basis which 
are generated from \aP) by applying to subsystem 2 the 
unitary transformation 



^2 1/3)2 ^e-'"^^ I 



(20) 



with ipp being one of the controllable phases of the in- 
terferometer. The resulting modified Werner states aw 
are obtained from Eq. (fT^ by replacing l^'^^^'') with the 
vectors 



|^lf)^(e-*^^|a/3)±e-«^°|/3« 



(21) 



while the corresponding modified twirling transformation 
is obtained by properly concatenating Twor of Eq. (jl7p 
with the unitary phase shifts transformations ([^0]) . i.e. 

fwerlp) = I dU {U® V2UV^) p (C/t ® V2U^V^) . 



In this context, the W{p) of Eq. (fT9|) is the maximum of 

(22) 



w{p)^Y.^¥-^\p\^^- 

a<f} 

for all choices of the phases {ipa}- 



B. Connection with the Fano factor 

In the notation of Eq. ^ the states \a(3) of Eq. ((12 
can be identified with 



(23) 



(this is possible since ai and a2 refer to independent an- 
nihilation operators). It is then easy to verify that, as- 
suming p to be the input state ^ and exploiting the nor- 
malization condition of the two-electron amplitude $a/3j 
the right hand side of Eq. (|22p becomes 



1 — w$ 



(24) 
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with being the quantity in Eq. ([8|). In particular, 
setting all the interferometer phases {(fa} to zero, we 
get W{p) of Eq. (fTB]) . Equation ((24|) has been explicitly 
derived for pure input states It can however be gen- 
eralized to any mixture by linearity. Combining Eqs. (|24[) 
and (O one gets 

This shows that by measuring the Fano factor one can 
determine and hence derive a lower bound for the 

entanglement of formation through Eq. (fT9|) . 



where |0) is now the electromagnetic vacuum (notice the 
absence of the productory with respect to E). Analo- 
gously, the current noise ([5]) is replaced by the function 



— / dh dt2 {Sh{ti)Sh{t2)), (27) 



which measures coincidence counts fluctuations of the 
photo-detectors located at the output ports 3 and AM 
Replacing Eq. ([26]) into Eq. jST]) gives 



V. CONCLUSION 

In this paper we introduced the fermionic analog of 
the Hong-Ou-Mandel interferometer. Within this con- 
text we discussed the possibility of detecting the presence 
of multi-mode entanglement for a class of two-electron 
inputs states by performing current noise measurements. 
Following Ref. 110 we also showed that Fano factors mea- 
surements provide a natural way of lower bounding the 
entanglement of formation of such incoming states. 

We end with a comparison of the results discussed 
above with the standard optical HOM interferometer— 
(see also Ref. [2l|) . The optical analog of the two-electron 
input state ([6|) is formally obtained by replacing aj,a{E) 
with bosonic annihilation operators aj{k) associated with 
coUinear electromagnetic modes of frequency Uk and 
propagating along the j-th optical pathj^ i.e. 

|*> = E *fci,fc.«I(fci)4(fc2) |0) , (26) 

ki ,k2 



F34 ^ -7== = -T(l - r)(l + w^) , (28) 

(-'3/\-'4) 

which should be compared with Eq. ^ of the electronic 
case (here w$ is as in Eq. ([8|). The sign difference in front 
of the function is a typical signature of the bunch- 
ing behavior of bosonic particles f2i Nevertheless also in 
the bosonic case one can use the negativity of as a 
signature of entanglement for the incoming two-photons 
states CT . 
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